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Abstract 

We present a phenomenological model of the quark-hadron transition in 
neutrino-nucleon scattering. Using recently extracted weak nucleon transition form 
factors, we investigate the extent to which local and global quark-hadron duality is 
applicable in the neutrino F±, F2 and F3 structure functions, and contrast this with 
duality in electron scattering. Our findings suggest that duality works relatively 
well for neutrino-nucleon scattering for the Fi and F3 structure functions, but not 
as well for F\. We also calculate the quasi-elastic, resonance and deep inelastic 
contributions to the Adler sum rule, and find it to be satisfied to within 10% for 
0.5 < Q 2 < 2 GeV 2 . 



1 Introduction 



Historically, neutrino scattering has provided vital information on the structure of the 
nucleon, complementary to that obtained by the more ubiquitous electromagnetic probes. 
In deep inelastic scattering (DIS), neutrino-induced structure functions have been used, 
in conjunction with electromagnetic structure functions, as the primary tool to separate 
valence and sea quark distributions. Neutrinos are also necessary to complete our knowl- 
edge of the full vector and axial vector structure of the nucleon elastic and transition form 
factors. 

At the parton level, deep inelastic structure functions describe incoherent scattering of 
a hard probe from quarks and gluons (generically, partons) in the nucleon; form factors, 
in contrast, characterize the coherent or bound-state response of the nucleon to an elec- 
tromagnetic or weak probe. While on the face of it the physics of coherent and incoherent 
processes is rather distinct, they are in fact intimately related through the phenomenon 
of quark-hadron duality. 

Quark-hadron duality in structure functions refers to the observation, first made by 
Bloom and Gilman PQ, that the average over resonances produced in inclusive eN scat- 
tering closely resembles the leading twist (or "scaling") function measured in the deep 
inelastic region. Furthermore, as Q 2 increases the average over resonances approaches 
the asymptotic scaling function. Within QCD, the degree to which this "Bloom-Gilman 
duality" holds is a direct reflection of the size of higher twist effects in the nucleon E] . 
According to the operator product expansion (OPE), higher twists are related to nucleon 
matrix elements of multi-quark or quark-gluon operators, which contain information on 
long-range, nonperturbative interactions between partons. Such interactions characterize 
the structure of the resonances and diminish with powers of 1/Q 2 as Q 2 — > oo. 

Recently there has been a resurgence of interest in duality in electron scattering at 
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Jefferson Lab and elsewhere, where its target, flavor, spin and nuclear dependence has 
been explored [H El El El El EUl HI] ■ Duality has been confirmed to good accuracy for the 
proton F2 and Fl structure functions to Q 2 values as low as 1 GeV 2 or even lower. The 
basic features of duality have also been studied in terms of dynamical models [121 EH HI] , 
in phenomenological parametrizations of the form factors ^Hl EH E] > as well as in the 
Rein-Sehgal model JH] and (in the A-resonance region) the Sato-Lee model j!9j . 

Within the models, neutrino scattering can provide an important consistency check, 
and lead to a better understanding of the systematics of nucleon iV — > resonance R 
transitions. While the phenomenological information on duality from electron scattering 
has been steadily accumulating ,20, , there is at present almost nothing known empirically 
about the workings of duality in neutrino scattering. There are plans, however, to measure 
neutrino cross sections using a high-intensity neutrino beam at Fermilab [2*Tj . 

In a parallel development, recent theoretical work has investigated the excitation of 
resonances by neutrinos for both J = 3/2 [221 EH] and J = 1/2 resonances [21]. In the 
latter work the weak vector form factors were determined from Jefferson Lab data using 
the conserved vector current (CVC) hypothesis, and two of the axial form factors from 
PCAC, although the Q 2 dependence was not very well constrained. To date there are 
only rudimentary data on neutrinoproduction of resonances beyond the P33(1232) region, 
however, more accurate data are expected, and a precise comparison will be possible 
in the future. In this paper we use the recent theoretical results to perform a detailed 
phenomenological study of duality in neutrino scattering. 

If one assumes that duality holds for neutrino scattering, then the average area under 
the resonances must follow the scaling curve. In this case the results of our comparison 
can be interpreted as a check on how well the Q 2 dependence of the transition amplitudes 
N — * R is known. Deviations from duality would in this case provide information on the 
size of the background, and of the axial form factors, which were not determined in the 
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model j2I] (for example, the normalization of C% and Cf for the _Di3(1520) resonance, as 
well as the Q 2 dependence of all axial form factors). Obtaining a better understanding 
of the dynamics in this kinematic region is also crucial for the interpretation of neutrino- 
oscillation experiments 

In Section |21 we review the formalism which is used in this study, and provide details 
about the transition form factor parameterizations. Results on local and global aspects 
of duality in neutrino scattering are discussed in Sec. EH and contrasted with duality 
in electron scattering. We also discuss the saturation of the Adler Sum rule, including 
its contributions from resonances, quasielastic and DIS regions. Finally, in Sec. 0] we 
summarize our results and draw conclusions from our study. 



2 Formalism 



Testing the degree to which duality in lepton-nucleon scattering is valid requires knowl- 
edge both of structure functions in the resonance region, and of the scaling functions 
applicable in the DIS regime. The former are calculated in terms of the nucleon — > res- 
onance transition form factors, while the latter can be evaluated from twist-two parton 
distribution functions. In this section we review both of these inputs, firstly outlining 
the parameterizations of the N — > R transition form factors from which the resonance 
structure functions are computed, and then summarizing the essential formulas for the 
twist-two structure functions. A more complete account of the formalism can be found in 
Refs. (2311211; here we shall present only those details which are pertinent to the specific 
discussion of duality. 
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2.1 Weak transition form factors 

In recent work by the Dortmund group, neutrinoproduction of the P33(1232) A resonance 
[221 123] was extended to cover also the second resonance region j21] , which includes three 
isospin-1/2 resonances: the Pn(1440) Roper resonance, and the two negative-parity states 
Di3(1520) and SS. In the following we summarize the weak transition form factors for 
these resonances. The definitions and notations for the cross sections and transition form 
factors are taken from Eqs. (IV.12)-(IV.15) and (IV.26)-(IV.28) in Ref. [23- 

P33(1232) resonance 

Historically, the P33(1232) (A) isobar has been studied more extensively than any other 
nucleon resonance. Electroproduction data on differential and integrated cross sections 
have been used to extract the N — > A transition form factor, and the resulting vector form 
factors, in the region Q 2 < 3.5 GeV 2 , can be parameterized (in the notation of Ref. [2*4] ) 

as 

w 2.13 D v {p) -1.51 D v w _ 0.48 D v 

3 ~l + Q 2 /4My 4 1 + Q 2 /AM 2 ' ° 5 " 1 + Q 2 /0.776M2 ' U 

where Dy = 1/ (1 + Q 2 /M V ) 2 is the dipole function with the vector mass parameter My = 
0.84 GeV, and the superscript (p) denotes a proton target. From isospin invariance, the 
electroproduction amplitudes of any isospin-3/2 resonance, are equivalent for proton 
and neutron targets, so that A(^n — > R^°) = A(^p — > R^ + ). Since the amplitudes 
are linear combinations of the form factors, the proton and neutron electromagnetic form 
factors are therefore also equal, C\ n ' = C\ p \i = 3,4,5. 

The weak vector form factors C\ (Q 2 ) for the amplitude A(W + n — > R^ + ) are related 
to the electromagnetic form factors. For isospin-3/2 resonances, these in fact coincide, 

Cf = Cf } = Cf , . = 3,4,5. (2) 
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The amplitude for a proton target is related to the neutron amplitude by Clebsch-Gordan 
coefficients, A(W + p -> R^ ++ ) = Vs A{W + n -> R^ + ). 

The form factors and Cj^ from Eq. (JTJ agree to within 5% with those obtained 
earlier under the assumption of magnetic dipole dominance. Since they are deduced from 
data for Q 2 < 3.5 GeV 2 , their normalization and Q 2 dependence should be reliable in this 
region. 



The axial form factors are obtained form PCAC, 

A _ 1.2 D A r A _ M 2 c 5 , o > 



where Da = 1/(1 + Q 2 /M\) 2 is the dipole term with the axial mass Ma = 1.05 GeV. For 
the other axial form factors, C^ 4 , we use the relations 

Ct{Q 2 ) = ~\c£{Q 2 ) and C£ = 0, (4) 

suggested by dispersion relations [213 HZ] • 

The P33(1232) resonance is known to be dominant for low energy neutrino scattering. 
The higher-mass resonances are very small for E u < 1.5 GeV and produce a noticeable 
peak in the invariant mass distribution for E u > 2 — 3 GeV. The second peak is produced 
primarily by the D13 and Sn resonances. 

-Di 3 (1520) resonance 

Among the isospin-1/2 resonances, RW, the L> 13 (1520) gives the dominant contribution 
in the second resonance region. The proton form factors in this case differ from those of 
neutrons. The vector part of the weak amplitude can be related to the electromagnetic 
amplitudes by isospin symmetry, 

A v (W + n -> R<V + ) = A^n -> R&+) - A( 1P -> R w+ ) . (5) 



Similarly, the weak vector form factors can be related to electromagnetic ones via 

Cr = C^-C^, . = 3,4,5. (6) 

The Q 2 dependence of the vector form factors (for Q 2 < 3.5 GeV 2 ) was determined in 
Ref. [21] from precise electromagnetic data from JLab in the second resonance region 

|2H1 123 ED], 

2.95 D v „m -1.05 Di 

V 



v 



23x3(1520): Cf = ^ " v 2 , C? = ^ 2 , <#> = -0.48 Z> 

6y ' 3 1 + Q 2 /8.9M 2 ' 4 1 + Q 2 /8.9M2 ' 5 

(w) _ -1.13 £> y (w) _ 0.46 £V r W__ ni7n 

° 3 "i + g 2 /8.9M 2 ' ° 4 "i + gv8.9M2' ° 5 " [) - UVy 



(7) 



for protons and neutrons, respectively. 



The normalization of the axial form factors is determined by PCAC and decay rates 
of the resonances. Unfortunately, their Q 2 dependence cannot be determined from the 
available data. In practice, we therefore consider two cases: (i) "fast fall-off", in which 
the Q 2 dependence is the same as for the P33 resonance, 



1 + Q 2 /3M 2 A ' 6 m 2 n + Q 2 



C i = 1 I ^2/0^2 > C t = M 2 ,^2 (" fast f all-off") ; (8) 



and (ii) "slow fall-off", in which the Q 2 dependence is flatter and has the same form as 
that for the vector form factors for each resonance, 

c ' = 1^ ' ci = M2 ^? ( " slow tall ' off " ) ' (9) 

The other two form factors, are unknown, and for simplicity we set them to zero, 
C 3 A = C£ = 0. 



Pxx(1440) and 5 , 11 (1535) resonances 

The two lowest-lying spin-1/2 resonances, Pn(1440) and 5*11(1535), both have isospin 
2 = 1/2. Their electromagnetic interaction depends only on two nonzero form factors, g% 
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and #2- For the proton these are determined for Q 2 < 3.5 GeV 2 from electroproduction 
helicity amplitudes, in analogy with the Di 3 resonance, 

2.3 Dy 



Pn(1440) : g\ 



(p) 



1 + Q 2 /4.3il4 ' 



and 



511(1535) : g\ 



(p) 

92 



(p) 



-0.76 D 



v 



1 - 2.8 In 1 + 



Q 2 



1 GeV^ 



2.0 D 



v 



1 + Q 2 /1.2M V 



1 + 7.2 In 1 



Q 2 



1 GeV^ 



9? 



0.84 D\ 



1 + 0.11 In 1 + 



Q 2 



(10) 



1 GeV 2 

For the neutron case, if one neglects the isoscalar contribution to the electromagnetic 



current, one can use the relation Wy 2 



-A^L In this case the general relation in Eq. 



between the weak and electromagnetic isovector form factors gives g. 



v 



-2ft 



(p) 



1,2. 



The axial vector form factors of these two resonances are constrained by PCAC, 



a a M(M R ± M) 

93 — 9i 



'12) 



Q 2 + ml 

with the ± corresponding to the Pu and Sn resonances, respectively. At Q 2 = the 
couplings are also determined from PCAC and the elastic vertex of the resonance decay, 
which is known from experiment, 

-0.51 D A 



P n (1440): gf(Q 2 ) 
9?(Q 2 ) 



1 + Q 2 /3M 2 



-0.51 D 



A 



Q 2 /4.3M 



;"fast fall-off") , 
;"slow fall-off") 



(13) 



and 

5 X1 (1535) : 9l (Q 2 ) 

9?(Q 2 ) 



-0.21 D A 

1 + Q2/3M 2 
-0.21 D A 
1 + Q 2 /1.2M 2 



"fast fall-off") , 
1 + 7.2 In ( 1 + 



Q 2 



1 GeV z 



"slow fall-off" 



(14) 
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2.2 Leading twist structure functions 



The second set of inputs required for duality studies are the inclusive structure functions 
F\ = MWi, F 2 = uW 2 , and F3 = vW% which describe the DIS region. Here we summa- 
rize the relevant expressions for the structure functions in terms of leading twist parton 
distribution functions (PDFs). In practice we use several PDF parametrizations, namely 
from the GRV [H], CTEQ 321 and MRST groups. 

For electron scattering, the F 2 structure function of the nucleon, defined as the average 
of the proton and neutron structure functions, is given (at leading order in a s and for 
three flavors), by 

Ff = \ (FT + Fn = %(u + u + d + d+ 2 -s + 2 -^j , (15) 

where the quark distributions are defined to be those in the proton. For neutrino scatter- 
ing, the corresponding F 2 structure function is given by 

F 2 N = x(u + u + d + d + s + s) . (16) 

In the moderate and large- 2 region, where strange quarks are suppressed, the weak and 
electromagnetic F 2 structure functions approximately satisfy the "5/18 rule", 

5x / -x 5 



F 



v - -(u + u + d + d) « — F% N . (17) 



The experimental confirmation of the factor 5/18 was indeed one of the important mile- 
stones in the acceptance of the description of DIS in terms of universal PDFs. 

In the Q 2 —>■ 00 limit, the Fi structure function is related to F 2 via the Callan-Gross 

relation, F 2 = 2xFi. Deviations from this relation arise due to perturbative a s corrections, 

as well as from target mass effects and higher twists. It is sometimes convenient also to 

define the longitudinal structure function F L , 

4MV, 
1 + ) 1 2 - 2./-K . 
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For large Q 2 the target mass term proportional to M 2 /Q 2 can be omitted; however, at 
Q 2 ~ few GeV 2 it can make an important contribution, especially at large x. Because 
the extraction of Fl requires longitudinal-transverse separation of cross section data, 
which is challenging experimentally, in practice the Fl structure function is not very well 
determined. For Fl we use the parametrization of the MRST group [33] • To estimate the 
uncertainty in its determination, we consider two different scenarios for 2xF±, namely (i) 
Callan-Gross relation, 2xFi = F2, and (ii) the exact expression for 2xF\ from Eq. p8j1 . 

Finally, the charge-conjugation odd F 3 structure function for neutrino-nucleon scat- 
tering is given by 

xF v z v = 2x(d-u + s), xFp = 2x(u - d + s) . (19) 

If one neglects the contribution of the strange quarks, s ~ 0, which is the case we con- 
sider here, the isoscalar Fg N structure function is given simply by the valence u v and d v 
distributions: 

xF% N « x{{u -u) + (d- d)) = x{u v + d v ) . (20) 

In the large-x region, where contribution of all sea quarks is very small, the F2 structure 
function will also be proportional to xF^, 

xF» N w F% N w ^Ff . (21) 

In the next section we will consider duality both for the total structure function, and 
for the valence-only structure function. In the context of "two-component duality" [37] , 
the resonance contributions are taken to be dual to valence quarks, while the nonresonant 
background is dual to the sea. In the resonance region, and especially at low Q 2 , it 
may be reasonable that a resonance-based model of structure functions would generate 
a valence-like scaling function. Indeed, there were strong suggestions of such resonance- 
valence duality in the recent proton F 2 data from JLab j3]. In the present study we test 

9 



the "two-component" duality hypothesis by comparing the calculated resonance structure 
functions with both the total and valence-only structure functions. 

In the next section we will use the above expressions to quantify the degree to which 
the averaged resonance structure function duals the leading twist structure functions for 
neutrino scattering, and compare this with duality for the electron case. 

3 Duality in electron & neutrino structure functions 
3.1 Electron scattering 

Before proceeding with the discussion of duality in neutrino scattering, we first consider 
duality for the better known case of electron scattering. Recent high-precision experiments 
at Jefferson Lab and elsewhere [H El El III El El EE UJ] have allowed accurate tests to be 
performed of Bloom-Gilman duality in electron scattering. For the proton F2 structure 
function, Niculescu et al. j3] found that the structure function in the resonance region, 
averaged over several intervals of x corresponding to the prominent resonance regions, 
reproduces well the scaling structure function down to relatively low values of Q 2 . Our 
aim here will not necessarily be to reproduce accurately the data with our resonance model 
[22 EH El] > but rather using phenomenological information on transition form factors to 
compare the workings of duality for neutrinos and electrons. 

The isoscalar nucleon structure function F 2 e7V = (F 2 ep + F 2 e ™)/2, calculated as a sum 
of electroproduced resonances, is displayed in Fig. ^ as a function of the Nachtmann 
variable £ = 2x/ yjl + 4M 2 x 2 /Q 2 for several values of Q 2 from 0.2 to 2 GeV 2 . The use of 
the Nachtmann variable takes into account kinematical target mass corrections, which can 
be important at large x and low Q 2 . The prominent peaks correspond to the P33(1232) 
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(A) resonance at the largest £ values in each spectrum. The next peaks, at smaller £, 
correspond to the second resonance region, where the 5*11(1535) and _Di3(1520) resonances 
dominate, and the Pn(1440) resonance gives a small contribution. With increasing Q 2 , 
the resonance peaks decrease in height and move to larger £. 

To examine the extent to which "local duality" works, we compare the £ and Q 2 
dependence of the individual resonances with the £ dependence of the leading twist F£ N 
structure function, for both the total and the valence-only cases. For the latter, we use 
leading twist PDFs at Q 2 = 10 GeV 2 from the GRV [HI], CTEQ 32, and MRST [33] 
groups. On average the resonances appear to oscillate around and slide down the leading 
twist function, reminiscent of the general features of the data as a function of £ and Q 2 
- see Refs. jUI^O]. For the calculated structure function, we consider the four resonances 
mentioned above, and integrate over the region 

1.1 <W< 1.6 GeV, (22) 

where the upper bound covers the range of the resonances taken into account in this 
analysis. 

The degree to which local duality is valid can be quantified by considering the ratio 
of integrals of the resonance (res) and leading twist (LT) structure functions, 



W) = J r - ' , - , (23) 



where Ti denotes F 2 , 1xF\ or xF 3 , and the integration limits correspond to £ m ; n = ^{W = 
1.6 GeV, Q 2 ) and £ max = £(W = 1.1 GeV, Q 2 ). The closer this ratio is to unity, the 
better the agreement with duality will be. Defining the ratio Ii{Q 2 ) in terms of integrals 
over the Nachtmann scaling variable £ instead of Bjorken x implicitly includes target 
mass corrections in the structure functions |3*H 1351 l3*o] . which are important at large x 
and small Q 2 . This is especially so for the Fl structure function, which is intrinsically 
small. An alternative approach would be to express the target mass corrected structure 
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% Q 2 , GeV 2 

Figure 1: Duality for the isoscalar nucleon F% N structure function. (Left) F% N as a 
function of £, for Q 2 = 0.2, 0.5, 1 and 2 GeV 2 (indicated on the spectra), compared with 
several leading twist parameterizations [3TJ IH21 EB] (valence and total) at Q 2 = 10 GeV 2 . 
(Right) Ratio I% N of the integrated F% N in the resonance region to the leading twist 
functions (valence and total). 

functions in terms of x and Q 2 [31] and perform the integrations over x [20] • For a first 
investigation of duality, and since we are mostly concerned about the relative differences 
between duality in neutrino and electron scattering, the integrals over £ in Eq. ()23j) will 
provide a sufficient test of integrated duality. 

The ratio 7| f° r electron scattering is shown in Fig. (right panel). The results 
are similar to those of the empirical analysis of JLab proton data 4J. The integrated 
resonance contribution is smaller than the leading twist at low Q 2 , but increases with 
increasing Q 2 . For Q 2 > 1 GeV 2 , the ratio If is within ~ 20% of unity when using the 
total DIS structure function. On the other hand, for the valence-only structure function 
the ratio is within ~ 20% of unity over a larger range, Q 2 > 0.5 GeV 2 . The better 
agreement of the resonance curve with the valence-only leading twist curve supports the 
notion of two-component duality [37 , as observed in the JLab _F 2 ep data |4J. In more 
refined treatments one would also take into account the Q 2 evolution of the leading twist 
structure function. This will modify the quantitative behavior of the ratio with respect 
to Q 2 , but not its essential features. 
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The fact that If < 1 in our model can be understood from the fact that only the 
first four resonances are included in the structure function. Since I2 is positive, the 
contribution from higher resonances as well as the nonresonant background will increase 
the numerator in the ratio If an d thus improve the accuracy of duality. The behavior 
of the ratio If N at large Q 2 is less well constrained due to the current poor knowledge of 
the leading twist structure function at high x and of the transition form factors at large 
Q 2 - 

Recently, new high-precision data from Jefferson Lab have allowed longitudinal- 
transverse separations to be performed, which have enabled the proton 2xF\ structure 
function to be accurately determined at large x jZj . This has made it possible for the first 
time to perform quantitative tests of duality for the F\ (or Fl) structure function. In 
Fig. |21 (left panel) we plot the isoscalar nucleon structure function 2xF£ N , calculated for 
the above mentioned four resonances, and compare with the leading twist parametrization 
from Ref. at Q 2 = 10 GeV 2 . The two leading twist curves correspond to the two sce- 
narios for 2xFi discussed in Sec. 12.21 namely, using the Callan-Gross relation, F 2 = 2xF 1 , 
and using the exact expression in Eq. (|18|). The difference between the two curves is 
relatively small at Q 2 = 10 GeV 2 , so that one can use either in the comparison with the 
resonance structure function. 

With increasing: Q 2 , the resonance 2xFf structure function is seen to slide along the 
leading twist curve, just as in the case of F£ N , but on average sits slightly higher than 
the leading twist curve. This can be quantified by considering the ratio If N , defined in 
Eq. (J23|) . which we plot in Fig. |2] (right panel). For most of the range of Q 2 > 0.5 GeV 2 
the ratio is some 30 — 50% above unity, which may indicate the need for additional terms 
in the resonance sum. On the other hand, it is known that target mass corrections have a 
relatively larger effect on 2xFi than on F2 |20~1 134j . and would tend to increase the leading 
twist functions, especially at large x (or £), and hence reduce the ratio If N . 
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MRST 
MRST(CG) 

.5 1 1.5 2 

Q 2 , GeV 2 

Figure 2: Duality for the isoscalar nucleon 2xFf N structure function. (Left) 2xFf N as 
a function of £, for Q 2 = 0.2,0.5,1 and 2 GeV 2 (indicated on the spectra), compared 
with the MRST parameterization ;33j at Q 2 = 10 GeV 2 , using Eq. (jl8|) (dotted) and the 
Callan-Gross (CG) relation, Fi = 2xFi (dot-dashed). (Right) Ratio If N of the integrated 
2xFf N in the resonance region to the leading twist function [HS] (see text). 

3.2 Neutrino scattering 

In the previous section we have demonstrated that the resonance model used here (23 
reproduces the qualitative features of duality observed in electron scattering, and have 
established the accuracy with which this duality holds in the model. Here we turn to 
the main aim of our paper, which is to compare and contrast the workings of duality in 
eiV scattering and in uN scattering. To make the comparison as rigorous as possible, 
we calculate the neutrino structure functions using the same four resonances as for the 
electron case, Figs. HH21 

Neutrino interactions have particular features which distinguish them from electro- 
magnetic probes. For the charge current reaction u^p — > A ++ , for example, only 
isospin-3/2 resonances are excited, and in particular the P33(1232) resonance. Because of 
isospin symmetry constraints, the neutrino-proton structure functions (F% p , 2xF" p and 
xF^ p ) for these resonances are three times larger than the neutrino-neutron structure 
functions. In this case the resonance structure functions are significantly larger than the 
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leading twist functions, F { > F { , and quark-hadron duality is clearly violated 
for a proton target. 

In neutrino-neutron scattering, in addition to isospin-3/2 resonances, all the isospin- 
1/2 resonances can also be excited. However, the total contribution of the three isospin- 
1/2 resonances considered here is smaller than that from the leading P33(1232) resonance. 
The leading twist curve for the vn structure functions lies above the resonance structure 
functions, F" n ( Tes > < _p^ n ( LT ) ; so ^hat quark-hadron duality does not hold for this case 
either. 

The general feature of the resonance curves is that at the onset of the resonance region, 
W < 1.6 GeV, the neutrino-proton structure functions are larger than the corresponding 
neutrino-neutron ones. In the deep inelastic region, on the other hand, the structure 
functions are larger for neutrino-neutron scattering. It has been argued ^2] for the case 
of electron scattering that for duality to appear one must sum over a complete set of even 
and odd parity resonances. In neutrino scattering, due to isospin symmetry constraints, 
duality will not hold locally for protons and neutrons separately, even if several resonances 
with both even and odd parities are taken into account [13] . In this case one can consider 
duality for the average of proton and neutron structure functions, which is the approach 
we take in this work. We will demonstrate below that in this case duality holds with even 
greater accuracy than for electron scattering. 

This discussion raises the question of how the transition occurs from the resonance 
to DIS regions in the case of neutrino scattering. We can speculate about the possible 
mechanisms of how this takes place. Starting from low W, the first resonance is the 
P33(1232), whose contribution to the up structure function is three times larger than 
that to the vn, as mentioned above. To compensate its influence, this resonance must be 
followed by several isospin- 1/2 resonances, which can only contribute to neutrino-neutron 
structure functions. This is what indeed happens — the Pn(1440), D 13 (1520) and SS 
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resonances are the next ones in the mass spectrum. In fact, the results of Ref. show 
that the P33(1232) form factors fall faster that a dipole, while those for the Z?i3(1520) 
and 5*11(1535) resonances fall slower than a dipole, so that at Q 2 « 2 GeV 2 the two peaks 
are comparable. From our calculations we also know that with only these resonances the 
vn cross section (and structure functions) are still smaller than those of vp. Additional 
resonances with higher masses may also follow this trend and further enhance the vn 
structure functions to overcome those for vp. 

At higher masses the isospin-3/2 resonances P33(1600) and 531(1620) appear, which 
also give three times larger contributions for vp scattering than for vn. They are again 
followed by the three isospin-1/2 resonances, Sn(1650), Di 5 (1675) and Fis(1680), two of 
which have spin 5/2. Their contributions can be large due to the summation over six 
final spin states, which further increases the neutrino-neutron structure functions. One 
could suppose that in this region the neutrino-neutron contribution would exceed the 
neutrino-proton. Furthermore, we have again one isospin-3/2 resonance, the .D33(1700), 
and three isospin-1/2 resonances — D 13 (1700), Pn(1710) and Pi 3 (1720) — to compensate 
its influence. Above W = 1750 MeV, and up to 2220 MeV, the isospin-3/2 resonances 
prevail, with 11 known states, and only 9 with isospin 1/2. 

A more detailed investigation of the interplay between the resonances with different 
spins would be possible after the form factors are determined for at least some of these 
higher- lying resonances. At present, however, we consider only the first four resonances, 
for which the vp cross section is always larger than the leading twist contribution, and 
the vn cross section is always smaller. This is one additional reason to compare only the 
average of the vp and vn structure functions. 

The neutrino-nucleon F% N structure function is displayed in Fig. |3] (left panel) as a 
function of £ for several values of Q 2 . Here the P33(1232) resonance is seen as the largest 
peak at each Q 2 . The next peak at lower £ (larger W) is dominated by the Di 3 (1520) and 
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SS resonances. The contribution from the latter becomes more significant with increasing 
Q 2 since its form factors fall off more slowly than the dipole. The contribution of the 
Pn(1440) resonance is too small to be seen as a separate peak. The two sets of resonance 
curves correspond to the "fast fall-off" (lower curves) and "slow fall-off" (upper curves) 
scenarios for the axial form factors discussed in Sec. 12.11 The smooth curves are obtained 
from Eq. (fT"""P using the GRV [3T] and CTEQ j32] leading twist parton distributions at 
Q 2 = 10 GeV 2 , as in Fig. [T] Just as in the case of electron-nucleon scattering, with 
increasing Q 2 the resonances slide along the leading twist curve, which is required by 
duality. As in Fig. 1, we show both the total structure function and the valence-only 
contribution. 
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Figure 3: Duality for the neutrino-nucleon F% structure function. (Left) F 2 i n the 
resonance region at several Q 2 values (indicated on the spectra), compared with leading 
twist parameterizations [3TJ E2] (valence and total) at Q 2 = 10 GeV 2 . (Right) Ratio 



tvN 



of the integrated _F 2 in the resonance region to the leading twist functions [3T*. 13*2] 
(valence and total). The upper (lower) resonance curves and the upper (lower) integrated 
ratios correspond to the "slow" ("fast") fall-off of the axial form factors. 



In Fig. [3] (right panel) we show the ratio of the integrals of the neutrino resonance and 
leading twist structure functions, defined in Eq. (|23)1 . The ratio is within ~ 20-25% of 
unity for Q 2 > 0.3 GeV 2 and, unlike the corresponding electron-nucleon ratio 7| JV does 
not grow appreciably with Q 2 . Again, the two sets of resonance curves correspond to 
the "fast fall-off" (lower) and "slow fall-off" (upper) scenarios for the axial form factors. 
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The difference between the curves reflects the uncertainty in the calculation of 1% • As 
expected, this ratio is close to 1 for the "valence-only" function at low Q 2 , which favors 
the hypothesis of two-component duality |37j . A comparison for Q 2 < 0.5 GeV 2 may 
be questionable, however, since there the perturbative QCD expansion is unlikely to be 
valid. For large Q 2 the ratio is sensitive to the parametrization used for the leading twist 
curve, and the difference between the two parametrizations is smaller than the difference 
between the valence-only and total functions. 

New features appear when considering the C-odd structure function Fg N . As discussed 
above for the case of F% , for the resonances considered here the proton F^ p structure 
function is larger than the neutron F% n , whereas for deep inelastic scattering the vn 
is larger. In our analysis we compare the isoscalar nucleon data, which are shown in 
Fig- IH (left panel). As before, the lower and upper curves in the second resonance region 
correspond to the "slow" and "fast" fall-offs of the axial form factors, respectively. 
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Figure 4: Duality for the neutrino-nucleon xFg N structure function. (Left) xF» N in the 
resonance region for several Q 2 values, (indicated on the spectra), compared with several 
leading twist parameterizations E2l E2] at Q 2 = 10 GeV 2 . (Right) Ratio I^ N of the 
integrated xF£ N in the resonance region to the leading twist functions. The upper (lower) 
resonance curves and the upper (lower) integrated ratios correspond to the "slow" ( "fast" ) 
fall-off of the axial form factors. 



To quantify the degree to which the resonance and deep inelastic structure functions 
are dual, we calculate the ratio of integrals for the xF^ N structure function as in Eq. (J23|) . 
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This ratio, shown in Fig. 0] (right panel), appears to fall off more rapidly with Q 2 than for 
the F£ N ratio, and reaches ~ 0.7 at Q 2 = 2 GeV 2 . The Fg N structure function is in general 
more sensitive to the choice of axial form factors, and our results are consistent with the 
uncertainty in the axial form factors, which is estimated to be ~ 30% at Q 2 = 2 GeV 2 . 




5 Q 2 , GeV 2 



Figure 5: Duality in the neutrino-nucleon 2xF" structure function. (Left) 2xF" N in the 
resonance region at several Q 2 values (indicated on the spectra), compared with the MRST 
parameterization (SSj at Q 2 = 10 GeV 2 using the exact expression in Eq.(jl8j) (dotted) 
and Callan-Gross relation (dot-dashed). (Right) Ratio 1% of the integrated 2xF" N in the 
resonance region to the leading twist function [33]. The upper (lower) resonance curves 
and the upper (lower) integrated ratios correspond to the "slow" ("fast") fall-off of the 
axial form factors. 

Finally, in Fig. El (left panel) we show the neutrino structure function 2xF^ N as a func- 
tion of £ for several Q 2 values. The resonance structure function are calculated for "slow 
fall-off" and "fast fall-off' axial form factors. The leading twist functions correspond to 
the MRST parametrization using the Callan-Gross relation and the exact expression 
in Eq.(|18|). As in the electron scattering case, the resonance contributions appear to lie 
above the leading twist curve for most of the range of £. The ratio I\ N , shown in Fig. El 
(right panel), is about 20% above 1 for Q 2 > 1 GeV 2 , which again may be an indication 
that target mass effects need to be removed from the leading twist structure function 
before comparing with the resonance contributions. 
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3.3 Adler sum rule 



One of the most fundamental results in neutrino scattering is the relation between the 
difference of the un and vp structure functions for quasi-elastic (QE) scattering and for 
the rest of the higher mass states PIU]. 

9iv E) (Q 2 )] 2 + [g?A E) (Q 2 )] 2 + [g% E \Q 2 )} 2 ^+J d » [ w Z n (Q 2 > ») - ^ P (Q 2 , u)]=2. 

(24) 

Because it measures the isospin of the target, this relation must hold for all values of Q 2 . 



In the Q — > limit, Eq. (|24|) is reduced to the Adler- Weisberger relation 
which has been verified experimentally to good accuracy. For Q 2 ^ 0, it is known as 
the Adler sum rule [40J, which has also been tested with data for neutrino deep inelastic 
scattering, and found to hold to ~ 20% accuracy 41]. At large Q 2 it has a simple 
interpretation in the parton model, in terms of integrals of valence quark distributions. 
Using the model [21] for the resonance form factors, we can study how the Adler sum rule 
is satisfied CIS db function of Q . 

For the QE form factors we use the following simple parametrization: 

9 r=-^-, s !r=^- m 

Dy Dy Da 

The W2 structure functions in Eq. ()24j) include contributions from resonance production 
and from the deep inelastic region. The resonance contribution is calculated for the 
first four resonances, as discussed earlier. The integration is performed in the range of 
^min < v < ^max corresponding to the final state mass range 1.1 < W < 1.6 GeV. In terms 
of £, the integration of the structure function for each Q 2 corresponds to the area under 
the resonance curve from £ min = £,{Q 2 ,W = 1.6 GeV) to £ max = £,{Q 2 ,W = 1.1 GeV). 
The contribution from the remaining £ interval, < £ < £ m in, corresponds to the higher 
W region. For this we assume that the structure functions are given by the leading twist 
contributions, calculated from the MRST parametrization [33] . 
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Figure 6: Decomposition of the Adler sum, as a function of Q 2 , into its QE (dashed), 
resonance (short dashed) and deep inelastic (dotted) contributions, as well as the total 
(solid). 

In Fig. |Ul the individual contributions from the QE, resonance, and DIS regions are 
plotted as a function of Q 2 . The (positive) QE contribution is large at low Q 2 but falls 
rapidly with increasing Q 2 . The resonant piece of the sum is negative, and partially 
cancels the QE component. The deep inelastic component grows with Q 2 , since £ — > 1 
as Q 2 — ► oo, and for Q 2 > 1 GeV 2 contributes some 80% of the integral. Combining 
the three terms, the sum rule is found to be satisfied within ~ 10% over the whole range 



Since the Adler some rule is based on very general grounds, one expects it to be exact. 
The 10% deviation of the calculated sum rule from the exact value should therefore be 
treated as an indication of the accuracy of the model. In practice, the uncertainty comes 
mainly from the axial form factors for the second resonance region, and suggests that 



0.5 < Q 2 < 2 GeV 2 . 
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some of them are underestimated. The requirement that the Adler sum rule is satisfied 
exactly could therefore serve as a restriction on the currently unknown axial form factors. 

4 Conclusion 

Motivated by the need to understand neutrino-nucleon cross sections in the Q 2 ~ few 
GeV 2 range, and the observation of quark-hadron duality in electron-nucleon scattering, 
we have performed a detailed phenomenological study of duality in neutrino structure 
functions. Using a recently developed model [221 I2H H3] for the first four lowest-lying 
nucleon resonances, we have computed the structure functions F2, 2xF\ and xF% in the 
resonance region for proton and neutron targets, and compared these with leading twist 
parameterizations. 

As a check of the resonance model, we have calculated the electron-nucleon structure 
functions and found that for each resonance these oscillate around the leading twist curves 
down to low values of Q 2 , in qualitative agreement with duality. For quantitative com- 
parisons, we defined ratios h{Q 2 ) of resonance to leading twist structure functions, which 
in the ideal case of duality, should be unity. Our results show that for the F% N structure 
function this ratio is below unity at low Q 2 , and slowly grows with Q 2 , consistent with 
recent experimental results The agreement with duality for 0.5 < Q 2 < 2 GeV 2 in 
this case is at the level of 20%. At low Q 2 the resonance averaged F 2 structure function 
resembles valence quark distributions, apparently oblivious to sea quark effects, which 
supports the hypothesis of two-component duality [37|. For the 2xFi structure function 
the ratio is about 40% above unity, but would be reduced after correcting for target mass 
effects in the leading twist structure function. 

For charged current neutrino scattering, duality does not hold for proton and neutron 
targets separately because of the dominant role played by the isospin-3/2 resonances. 
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However, averaging over proton and neutron targets leads to large cancellations between 
7 = 3/2 and 7=1/2 resonances, so that duality holds at the 20% level for isoscalar 
vN structure functions. Furthermore, the ratios I 2 N (Q 2 ) and I£ N (Q 2 ) appear to reach 
constant values already for Q 2 ps 1 GeV 2 . 

Another interesting feature of our analysis is that the ratios If N (Q 2 ) and I" N (Q 2 ) 
of the 2xF 1 structure functions are consistently above unity. This may be an indication 
of the importance of target mass corrections in the leading twist Fi structure functions, 
which are known to be more important than those in F 2 . 

In these comparisons we have used leading twist structure functions obtained from 
global parton distributions, which are well constrained by experimental data, especially 
for F 2 . For the resonances, on the other hand, the data are very sparse, and theoretical 
input needs to be used. Our results therefore have an inherent uncertainty arising from 
poor knowledge of the transition form factors, particularly at high Q 2 . 

The results obtained here raise the following question: what is the most efficient and 
quantitative method for comparing the resonance contributions with the scaling curves 
and their QCD corrections? One approach is to compare the various contributions to sum 
rules, in which integrals over resonance and DIS contributions must reproduce physical 
constants. To this end we computed the various contributions to the Adler sum rule as 
a function of Q 2 . This exercise shows how the relative contributions vary with Q 2 , and 
saturate ~ 90% of the sum rule. The remaining 10% could be accounted for by including 
more resonances, and by better determining the transition form factors. 

Overall, our quantitative study of neutrino reactions indicates that duality in structure 
functions, averaged over protons and neutrons, is expected to work to even better accuracy 
for neutrino scattering than for electron scattering. 
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